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1. Calcul tensoriel 


Exercice A. Montrer que la symetrie est une propriete tensorielle, c'est-a-dire que si un 
tenseur S est symetrique (S^ — Sj L ) dans une base orthonormee donne B — (eq, e^, e^), 
alors cette propriete est egalement vraie dans toute autre base orthonormee B' — 
(e'i, e' 2 /e 3 ). 

Exercice B. Montrer que l'antisymetrie est egalement une propriete tensorielle. 

Exercice C. Montrer qu'un tenseur T quelconque peut toujours se decomposer en une 
partie symetrique et une partie antisymetrique. 

Exercice D. Soit S un tenseur symetrique et A un tenseur antisymetrique, montrer que 
Ton a toujours 5: A — 0. 


Exercice E. Soit S un tenseur symetrique et T un tenseur quelconque, montrer que l'on 
a toujours : S: T — S\ T s , ou T s est la partie symetrique de T. 

Exercice F. Soit A un champ vectoriel. Montrer que l'on a toujours div(roi v4) = 0 

Exercice G. Soit A un champ scalaire. Montrer que l'on a toujours rot (grad A) = 0 

Exercice H. Soit la base curviligne polaire (e^ - , e^, e^), telle que : 


x — R ■ cosO ■ sincp 
y = R ■ sinO ■ sincp 
z — R ■ coscp 




1. Calculer la surface exterieure d'une sphere de rayon R, sachant qu'un element de 
surface infinitesimal peut s'ecrire dS — R 1 2 sincpd6dcp. 

2. Calculer l'integrale du champ coscp ■ sur cette sphere, sachant que le vecteur 
elementaire dans la direction radiale vaut = sincpcosde% + sincpsindey + coscp~e^ 
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Corriges 


Exercice A. 

Soit un tenseur S symetrique. Dans la base B — (e^e^e^) on peut done ecrire 
Sij — Sji . Soit une base orthonormee quelconque B r — {e' lt e' 2 >' e -$) differente de B. 
On appelle a la matrice de passage de B a B'. Si on appelle S\j les termes de la 
matrice de S dans la base B', on a d’apres le cours (en notation d'Einstein) : 

S'-■ = a- a- S 

*-> ij ^ip^jq^pq 

Or la matrice S est symetrique, on a done S pq — S qp . On ecrit done : 

S'-■ = a- a- S 

*“> ij ^ip^jq^qp 

Les indices p et q du second membre sont muets, on pourrait done les remplacer par 
n'importe quelle lettre, et on peut aussi les intervertir : 

S ij OCiqOCjpSpq ^jp^iq^pq 

On en deduit que S'\j = S'j h ce qu'il fallait demontrer. 

Exercice B. 

Soit un tenseur S antisymetrique. Dans la base B — (e^e^e^) on peut done ecrire 
S t j — —Sji. Soit une base orthonormee quelconque B' — (e' 1 ,e' 2 ,'e 3 > j differente de B. 
On appelle a la matrice de passage de B a B'. Si on appelle S\j les termes de la 
matrice de S dans la base B', on a d'apres le cours (en notation d'Einstein) : 

S'-■ = a- a- S 

*-> ij ^ip^jq^pq 

Or la matrice S est symetrique, on a done S pq — —S qp . On ecrit done : 

S' - = —a- a- S 

*“* ij Ll ip' A 'jq Lj qp 

Les indices p et qdu second membre sont muets, on pourrait done les remplacer par 
n'importe quelle lettre, et on peut aussi les intervertir : 

S ij ^iqOCjpSpq ^jp^iq^pq 

On en deduit que S'ij = — S'ji, ce qu'il fallait demontrer. 
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Exercice C. 

On reprend les formules du cours. Soient les tenseurs T s et T A donnes par les formules 
suivantes : 


t s = i(r + r r ) 

f A = ^(T-T T ) 

On va demontrer que T s est symetrique et que T A est antisymetrique. Pour cela on se 
place dans une base B — (tq, e^, e^), pour laquelle le tenseur T s'exprime sous la forme 
d'une matrice de terme generale T L j. Explicitions les termes des matrices de T s et T A : 

T l S ,=l(T„ + T Jl ) = ^(T, l + T l ,) = T? l 

T,1 = \ fa, - 7},) = - \ (T n - T tj ) = —Tji 
On a done demontre que T s est symetrique et T A antisymetrique. 


Exercice D. 

Dans une base orthonormee B — (e^, e^, eyT), on peut ecrire que S L j = et egalement 
que Aij — —Aji . On en deduit que, comme pour tout tenseur antisymetrique, les termes 

diagonaux de la matrice de A sont nuls dans toute base. D'apres le cours, le produit 
doublement contracts de ces deux tenseurs est un scalaire egal a : 

S:A=S tJ A Jt 

Developpons cette notation d'Einstein sous forme explicite : 


S'. A — + S-12A21 + S 13 A 31 

+ S 2 iA 12 + S 22 A 22 + ‘^ 23^32 
+ ^ 31^13 + S22A23 + ^ 33-^33 

On sait que l'on a A 1± = A 22 — A 33 — 0, et par ailleurs : 

“-*21 = S 12 ; S 31 = S 13 ; S32 - S 2 3 
A21 — —A 12 ; A 31 - —A 13 ; A 32 — —A 23 


On en deduit directement S : A — 0. 
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Exercice E. 

Le resultat decoule directement de celui de l'exercice precedent et de la distributivite 
de l'operateur produit doublement contracts : 

S- T = S: (f s + f ■") = S: T s + S: T* = S: T s 


Exercice F. 


Posons le resultat intermediaire B = rot A. Dans une base orthonormee donnee 
B — (el, el, e^), on a par definition : 




B = 


dA 3 

dA 2 

dx 2 

9x 3 

9A ± 

9A 3 

dx 3 

9x ± 

dA 2 

9A ± 

,dx x 

dx 2 


\ 


Par ailleurs, on a par definition : div B — — + — + —. 

r dx t dx 2 dx 3 


On en deduit que div rot A vaut: 

d (9A 3 dA 


div rot A — 


2 ' , J?_ (^1 _ + ( <>A 2 _ d£i\ 

Vdxo dxi) + dx\dxi dx 7 ) 


dA 2 9A ± ' 


dx 1 \dx 2 dx 3 J dx 2 \dx 3 dx 1 ) dx 3 \dx 1 dx 2 


div rot A = 


d 2 A, 


d 2 A 


2 + d A t 


d 2 A, 


+ 


d 2 A-: 


d 2 A 1 


dx t dx 2 dx 1 dx 3 dx 2 dx 3 dx 3 dx ± dx 3 dx 1 dx 3 dx 2 

L'ordre des derivations partielles successives d'une fonction de plusieurs variables est 

_^ —> 

quelconque, on peut done en deduire directement: div rot A — 0. 


Exercice G. 

On se place egalement dans une base orthonormee B — (el, el, el). Soit le resultat 
intermediaire B = grad A donne par la formule du cours : 

/—\ 

dx ± 

dA 

B ~ 

dA 

\dxj B 
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Par ailleurs, on peut ecrire, toujours dans la base B : 

/<3Z?3 dB 2 


rot B — 


On en deduit: 


\ 


rot grad A = 


dx 2 

dB ± 

dx 3 

dB 3 

dx 3 

dB 2 

dx x 

dB ± 

oT 

X 

I-* 

dxj 

( d dA 

d 

dx 2 dx 3 
d dA 

1 

QjI 

dx 3 dx ± 
d dA 

1 1 
1 

QJ X \ 


dA 


)x x dx 2 


\ 


3 

dA 


1 u *3 

dA 


Pour la meme raison que dans l'exercice precedent, on peut done ecrire directement la 
formule classique rot grad A = 0. 


Exercice H. 

1. Un element de surface dS de la sphere de rayon R delimite par deux secteurs 
d'angles infinitesimaux dO et dtp s'ecrit: 

dS — R 2 sirupd6d(p 

On cherche a calculer S = f ., dS. En parametrant la surface en fonction de Q et cp, 

J sphere 1 ” 

on peut expliciter cette integrale : 


rTC r2.1T 

5=1 | R 2 sincpddd(p 

J(p =o Je=o 

Le rayon R est independant de 6 et cp, et le terme sincp est independant de G, on peut 
done ecrire : 


5 = 



dtp 


S — R 2 f sin(p[6]l n dtp — 

J(p = o 



2tiR 2 [—cos(p]q 
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Finalement, on obtient le resultat classique S = 4nR 2 . 


2. On doit calculer : 


r r TC r 2.71 

A — I coscp-e^dS — I I coscp -e^ ■ R 2 sincpd9dcp 

Jsphere J<p =0 •'0 = 0 

Avec : = sincpcosOe% + sincpsinOey + coscpe1 

On est en coordonnees curvilignes, done on ne peut pas sortir le vecteur car il n'est 
pas independant du point d'integration (du point de la sphere de surface dS). En 
revanche, les trois vecteurs de la base cartesienne ont cette propriete, et peuvent etre 
sortis de l'integrale. On peut done remplacer ep par son expression, et ecrire dans la 
base B — (e^, e 2 , e^) : 


A = 


/ r 71 r2ll 

I J<p=oJg=o 

rTZ r2lZ 

j(p= 0 ■'0 = 0 

0 2 n 

) = 0 


coscp ■ sincpcosO ■ R 2 sincpd6dcp 


\ 


coscp ■ sincpsinO ■ R 2 sincpd6dcp 


\ 


coscp ■ coscp ■ R 2 sincpd6dcp 


) 


I I I cos(pcos0sin 2 (pd0d(p\ I I coscpsirdcp I I 

I J(p= 0^0=0 1 I '(^=0 \Je=o 

n r 2n 


A = R 2 


r 11 r Z 11 

J(D = 0 Jq = 0 


coscpsindsin 2 cpd0dcp 


rTC r 2.71 I I rTC 

\ I I cos 2 (psimpddd(p / \ I 

\ J(p=oJe=o / B \ J<p =0 

Du fait de la periodicite des fonctions trigonometriques, on a : 

/-27T r 2.71 

I cosOdO — I si 
'e =0 h =0 


= R 2 


coscpsin^cp l 1 cosOdO 
<p =0 \^0 =0 / 

n 


B 

2n 


jdcp\ 

I coscpsin 2 (p I I sinOcpdd ) dep 
lm = [) \jfl = n ) 


'0 = 0 

cos 2 cpsinq){ I d0)dcp 

(p =0 \”'0 = O 


/ 


sinddd = 0 


Par ailleurs, on a : 


r2n 

de = 

Je =0 


= 271 
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2. Cinematique et deformations 


Probleme. On considere un mouvement defini dans la base B — (e^,e^,e^) par sa 
representation lagrangienne (oj est une constante positive) : 

'x 1 — X^osicot) — X 2 sin(a)t ) 

■ x 2 — X^inicot) + X 2 cos(a)t ) 

, *3 - ^3 

1. Calculer le tenseur gradient F, le tenseur des dilatations C, et le tenseur des 
deformations E de ce mouvement au point X et a l'instant t. 

2. A quelle classe particuliere ce mouvement appartient-il ? 

3. Pour un instant t donne, calculer la dilatation en un point X et dans une direction dX 

—> 

4. Pour un instant t donne, calculer le glissement en un point X et pour deux directions 
orthogonales dX et dX’. 

5. On considere un milieu anime de ce mouvement, muni d'une masse volumique 
homogene p 0 a l'instant t 0 — 0. Calculer le jacobien de la transformation, ainsi que la 
masse volumique du milieu a l'instant t. 

6. Calculer le champ de vitesse V(X, t) et le champ d'acceleration y(X, t) en 
coordonnees lagrangiennes. 

7. Exprimer les coordonnees initiales a partir des coordonnees actuelles. Calculer le 
champ de vitesse V(x, t) et le champ d'acceleration y(x, t) en coordonnees 
euleriennes. 

8. Calculer les tenseurs des taux de deformations euleriens D(x,t ) et des taux de 
rotation £l(x, t). 

9. On definit les coordonnees polaires lagrangiennes ( R,Q,X 3 ) par le changement de 
variables (X 1 ,X 2 ,X 3 ) = (R ■ cos®, R ■ sin®,X 3 ) et les coordonnees euleriennes 
(r,0,x 3 ) par le changement de variables (x 1 ,x 2 ,x 3 ) — (r ■ cos0,r ■ sin0,x 3 ). 
Expliciter les fonctions X r et X g definissant une nouvelle representation lagrangienne 
du mouvement de la forme : 


(r — X r (R,Q,X 3 , t) 
19 =X e (R,Q,X 3 ,t) 
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10. On definit une base curviligne polaire B' — ( e(, e^ - , e^), telle que : 

( ep = ~elcos9 + ~e^sin6 
le q = —elsinO + ~e(cos9 

On note alors le champ de vitesse dans cette nouvelle base : 


V(x, t) = V r (r, 9, x 3 , t)T r + V e (r, 9, x 3 , t)e£ + E 3 (r, 9, x 3 , t)e( 

Indiquer l'expression des composantes polaires V r , V g et V 3 du champ de vitesse 
eulerien. 

11. Calculer l'acceleration centrifuge y r (r, 9, t ) et l'acceleration tangentielle y g (r, 9, t ) 
du mouvement etudie. 

12. Definir les trajectoires associees a ce mouvement 

Examen partiel : Etude cinematique d'un tourbillon. On considere un mouvement, 
defini dans la base orthonormee B — par la representation eulerienne 

suivante : 

'V\ = -A(x 1 ,x 2 )x 2 
■ V 2 — A(x 1 ,x 2 )x 1 

{ ^3 = 0 

Dans cette representation, A est une fonction scalaire des deux coordonnees x ± et x 2 , 
definie par l'expression A(x ± ,x 2 ) — ~(^i + ^|) _1 °u a est une constante positive. 
Ce mouvement est done uniquement defini pour || x\ + x\ [| > 0. 

1. De quel type de mouvement s'agit-il ? 

2. Montrer que la matrice du tenseur gradient de vitesse L = grad V dans la base 
B — (el, el, e() s'exprime par : 


= _ 

xi + x 2 


■ 2x x x 2 

2 2 
x 2 — x ± 

I- 0 


2 2 

x 2 — X-y 

-2XyX 2 

0 


O' 

0 

0- 


B 


3. En deduire immediatement les matrices du tenseur des taux de rotation fl(x, t ) et du 
tenseur des taux de deformation euleriens D(x, t ). 


4. On considere un milieu continu de masse volumique p 0 a l'etat initial (t = 0) et 
anime de ce mouvement. Calculer la divergence du champ de vitesse, et en deduire la 
masse volumique en tout point et a tout instant t > 0. 
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5. Montrer que l'acceleration y(x 1 ,x 2 ,x 3 ,t) en coordonnees euleriennes cartesiennes 
(c'est-a-dire dans la base B — (fq, e^, e^)) s'exprime par : 


y(x, t ) = A 2 (x 1 ,x 2 ) 


r~x ± - 

-x 2 

0 


B 


6. On definit les coordonnees polaires lagrangiennes (R,Q,X 3 ) par le changement de 
variables (X lt X 2 ,X 3 ) — (R ■ cosQ, R ■ sin®, X 3 ) et les coordonnees euleriennes 
(r, 9, x 3 ) par le changement de variables (x 1; x 2 , x 3 ) = (r ■ cosO, r ■ sinQ, x 3 ). 


On definit egalement une base curviligne polaire B' — (ep, e^, e^), telle que : 

( ep = e^cosO + e^sinO 
le g — —elsin9 + ~e^cos9 


Calculer l'acceleration y(r, 6, x 3 , t ) en coordonnees euleriennes polaires. 


7. Exprimer le champ de vitesse en coordonnees euleriennes polaires. En deduire que 
les particules ont des trajectoires circulaires autour de l'origine, de vitesse angulaire : 

a 


Donner sans calcul l'expression de la representation polaire lagrangienne du 
mouvement sous la forme r = X r (R, 0, x 3 , t) et 6 — X e ( R, 0, x 3 , t). 


8. Demontrer que la representation lagrangienne du mouvement dans la base 
cartesienne sous la forme x — x^X, t) s'exprime par : 


x = 


cos(ojt) —sin(ojt) 0 
sin(ojt) cos(ojt) 0 
0 0 1 




Montrer que, dans cette expression, on a : co - ^ = A(x 1 ,x 2 ). 


9. Definir les trajectoires associees a ce mouvement, ainsi que les champs de vitesse et 
d'acceleration 


10. On considere le point de vecteur position x = r ■ e^. Calculer la matrice du tenseur 
des taux de deformations euleriens en ce point. De quel type de deformation s'agit-il ? 
En deduire les valeurs propres et les vecteurs propres du tenseur des taux de 
deformations euleriens en ce point. 
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Corriges 


Probleme. 

d x ■ 

1. Le terme general de la matrice du tenseur gradient dans la base B est done 

on a : 


F = 


COSOt 

sinot 
. 0 


—sinot 0 
cosot 0 


Pour le tenseur des dilatations, on a = F pi F pj - en notation d'Einstein, done : 



cos 2 ot + sin 2 ot 
cosot ■ sinot — cosot ■ sinot 
0 


cosot ■ sinot — cosot ■ sinot 
cos 2 ot + sin 2 ot 
0 


0 " 

0 = / 
1-^B 


Enfin, on sait que E — ^ [C — 7 ), et on a done E — 0. 


2. Le tenseur des deformations est nul, on est done en presence d'un mouvement 
rigidifiant. 


3. Puisque C = I, on peut dire que toute direction est direction principale. La dilatation 
dans une direction quelconque dX vaut done A^dX^) — = 1 


4. Pour la meme raison, le glissement entre deux directions orthogonales quelconques 


dX et dX' vaut: y (dX, dX'^\ — , Clj =0 

V / IcuCjj 


5. Le jacobien de la transformation est le determinant de F. On a done : 


J — detF — 


cosot —sinot 0 
sinot cosot 0 
0 0 1 


= cos 2 ot + sin 2 ot — 1 


Par consequent la masse volumique du milieu est constante dans le temps et en tout 
point. 
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6. Le champ de vitesse V(X, t) s'obtient par derivation materielle du champ de 


vecteurs exprimant les positions actuelles des particules : V ( X, t) = P ^'^ . On a : 


r x 1 — X ± cos(cot) — X 2 sin(cot ) 
x 2 — X^inCcot) + X 2 cos(cot ) 

Xo — Xo 


Done : 


Dx ± , . 

V t = = co ■ (—X^inicot) — X 2 cos(cot)) 

Dx 2 , . 

V 2 — = at ■ (X^osicot) — X 2 sin(cot)) 


Dt 


Dxo 

V 3 — ——— — 0 
3 Dt 


De meme, le champ d'acceleration s'obtient par y(X, t) = DV ^ C ' t \ done : 


y 1 — = a) 2 ■ ^—X^ os (cot) + X 2 sin(cot )) 


<Y 2 = 


Dt 

PV 2 

Dt 


— co 2 ■ (— X^inicot) — X 2 cos(cot )) 


DV , 

y 3 = —- = 0 
n Dt 


Ces champs s'expriment en fonction de X , et sont done bien en coordonnees 
lagrangiennes. 


7. D'apres l'enonce, les coordonnees actuelles s'obtiennent a partir des coordonnees 
initiales par le systeme suivant: 


'Xi 


cosoot —sincot 0" 


r*ii 

*2 

— 

sincot coscot 0 


*2 

X 3 . 


o 

o 

h- 1 




On isole en particulier les deux premieres equations de ce systeme : 

r*il _ Icoscot —sinojtl \X ± 

\-X 2 \ \-sincot coscot J [X 2 

On cherche a inverser ce systeme pour exprimer X en fonction de x. On utilise la 
formule d'inversion d'une matrice 2*2 : 

la bV 1 _ 1 r d -hi 

I c di ad — be L—c a J 
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Done : 


rcosojt 

—sincotA 

1 _ 1 r coscot 

sincot 1 

[coscot 

—sincot 1 

Uinojt 

coscot l 

cos 2 cot + sin 2 cot sincot 

coscoti 

[sincot 

coscot J 


Cette matrice est orthogonale, car son inverse est egale a sa transposee. Elle traduit 
done une rotation. On en deduit: 


" X ± — x ± cos(cot) + x 2 sin(cot ) 

• X 2 = —x^inicot) + x 2 cos(cot ) 

, ^3 - *3 

Dans l'expression du champ de vitesse et du champ d'acceleration calcules a la 
question precedente, on peut alors remplacer les coordonnees initiates X par leur 
expression en fonction des coordonnees actuelles x. On obtient la formulation 
eulerienne suivante : 


{ V 1 — co ■ [x^os^cot) + x 2 sin(cot))sin(cot) — ( 

V 2 =a>- ((^ cos{cot ) + x 2 sin{cot))cos{co€) — (- 

E 3 = 0 


-x^in^cot) + x 2 cos(cot))cos(cot)^j 
x^in^cot) + x 2 cos{cot))sin{cot )) 


(Vi = -ux 2 
j V 2 = mx 1 
V E 3 = 0 



y 2 = co 2 


(^—(x 1 cos(cot) + x 2 sin{cot))cos{co€) + ( 

(^cosCuit) + x 2 sin(cot))sin(cot ) — (■ 

Ks = 0 


■x^inicot) + x 2 cos 
x^inCcot) + x 2 cos 


(cot))cos(cot)^j 


f 2 

y 1 — —co ■ x x 
■ y 2 = -co 2 ■ x 2 
, 73 = 0 


Ce dernier resultat pouvait aussi s'obtenir en passant par la formule donnant 
directement l'acceleration par y — ^ + (grad V^V, avec ^ = 0 et: 


grad V 


0 

co 


LO 


—co 0 
0 0 
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8. On calcule d'abord le tenseur gradient de vitesse L = grad V, avec la formule de 
cours (grad V).. = — L . On remarque qu'il s'agit d'un gradient eulerien (g minuscule), 

IJ OX j 

done les derivations sont effectuees par rapport aux coordonnees actuelles x : 


L — grad V — 


0 

CO 

LO 


—co 0 
0 0 


Ce tenseur est antisymetrique, on peut done directement ecrire D — 0 et H = L. Ces 
resultats peuvent etre retrouves par le calcul avec les formules de cours : 

D = i(L + l7) et n = i(L-I r ) 


9. En coordonnees polaires, on a les correspondances suivantes : 

(X 1 — R • cos® Cx 1 - r ■ cosd (r = X r (R, 0, X 3 ,t) 

{X 2 = R-sin® ’ [x 2 = r ■ sinQ ’ [6 = X g (R,Q,X 3 , t) 

On cherche les expressions des fonctions X r et X e . On peut ecrire a partir de ces 
expressions : 


r = 

A 

x 2 + x\ 

cosd — — 
r 

) 

sind — 

a 2 

r 

R = ) 

Xl + x\ , 

Xi 

cos 0 = — 

R 

t 

sin 0 = 

x 2 

R 


Or on sait d'apres renonce que, dans le repere cartesien, le mouvement s'exprime par : 

= X^os^cot) — X 2 sin(a)t ) 

■ x 2 — X^inicot) + X 2 cos(a)t ) 

, *3 - ^3 

On peut remplacer x ± et x 2 par leurs expressions dans la formules donnant r : 


x, 2 + x, 2 


r — 



On en deduit l'expression de la fonction X r . On effectue la meme operation avec 
l'expression de cosQ : 
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X t cos{oot) — X 2 sin(oot) 

cosQ = - 

r 

Dans cette derniere expression, on remplace X ± , X 2 et r par leur expression polaire : 

R ■ cosQcosioot ) — R ■ sinQsin(oot) 

cosQ — ---= cos(0 + go t) 

R 

On en deduit 0 = 0 + oot, ce qui nous donne la forme de la fonction X g . Finalement, 
on a: 

f r — R 
10 = 0 + oot 


10. On travaille maintenant dans la base curviligne polaire B' = (ep, ~eg, e^) telle que : 

( ep = e^cosO + 'e^sinO 
le q = —elsinO + ~e^cos6 

On cherche a exprimer le champ de vitesse eulerien dans cette base. On a demontre 
dans la question 7 que le champ de vitesse pouvait s'exprimer par : 


V — —oox 2 ■ e x + oox t ■ e 2 

Dans cette expression, on remplace x ± et x 2 par leurs expressions en coordonnees 
polaires : 

V — —o)(r ■ sinQ ) ■ ~e[ + oo(r ■ cos0) ■ e^ 


done : 


V — oor(—sinQ ■ ~e{ + cosQ ■ e£) — oor~e^ 

Dans la base B' — (ep, eg, e^), la vitesse s'exprime done par une formule tres simple : 



11. Dans le meme ordre d'idee, on cherche a exprimer la formulation eulerienne du 
champ d'acceleration dans la base B'. On a montre que le champ d'acceleration 
s'exprimait par : 

f = —oo 2 x 1 ■ ~el — oo 2 x 2 ■ e^ 

On remplace x x et x 2 par leurs expressions en coordonnees polaires : 
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f — —to 2 (r ■ cos0) ■ ~e{ — to 2 (r ■ sm0) ■ 
y = —n> 2 r(cos0 ■ + sin0 ■ e^) -- —to 2 r~e^ 

Dans la base B r — l'acceleration est done purement centripete : 



12. Les trajectoires des particules materielles sont des cercles centres sur l'origine. La 
vitesse angulaire est constante, la vitesse est purement orthoradiale (perpendiculaire a 
un rayon) et augmente proportionnellement a la distance a l'origine. L'acceleration est 
purement centripete (dirigee vers l'origine) et est egalement proportionnelle en norme 
a la distance a l'origine. 


Examen Partiel. 

1.11 s'agit d'un mouvement plan, car la vitesse dans la direction est nulle et les deux 
autre composantes V x et V 2 sont invariantes par translation dans la direction e^. Par 
ailleurs la representation eulerienne du champ de vitesse est constante dans le temps, 
done on est par definition en presence d'un mouvement permanent. 

2. le tenseur gradient de vitesse L = grad V s'exprime dans la base B — (e^, e^, e^) 

d V' 

par la matrice de terme general L £; - = On travaille d'abord sur la derivation de la 
fonction A(x 1 ,x 2 ) = — (x 2 + x|) _1 : 


dA(x ± ,x 2 ) 

a 

2x ± 

-Aix^xy ■ 

2x t 

dx ± 

2n 

(x 2 + x 2 ) 2 

x 2 + x 2 

dy4(x 1; x 2 ) 

a 

2x 2 

-Aix^xy ■ 

2x 2 

dx 2 

2n 

(x 2 + x 2 y ~ 

X 2 + X 2 


On peut alors calculer les termes non triviaux de la matrice de L, a commencer par les 
deux premiers termes diagonaux : 


dV 1 3 A(x 1 ,x 2 ) 

dx 1 2 dx t 


A(x t ,x 2 ) 


2x t x 2 
x 2 + xl 
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dV 2 _ dA(x 1 ,x 2 ) 

dx 2 1 dx 2 


-A(x lf x 2 ) 


2*1*2 

*1 + *1 


On a egalement deux termes non-diagonaux assez complexes : 


dV x 

dx 2 


~x 2 


8A(x 1 ,x 2 ) 

dx 2 


A(x i,* 2 ) 


A(x i,* 2 ) 


2 *| 

*1 + *f 


^(*i,* 2 ) 


dFj 

d* 2 


21(*i,* 2 ) 


2 2 
*2 — *1 

*1 + *2 


<^2 

d*i 


*1 ■ 


dy4(* 1; * 2 ) 

d* x 


+ j4(*i,* 2 ) 


-4(* lf * 2 ) 


2*i 

*1 + *2 


+ j4(*i,* 2 ) 


<^2 

d*! 


v4(*i, * 2 ) 


2 2 
*2 — *1 

*1 + *2 


Tous les autres termes de la matrice de L sont nuls car on est en mouvement plan. 
Finalement, le tenseur gradient de vitesse s'exprime sous forme matricielle dans la 
base B — (e^, e 2 , e^) par : 


Y _ ^(*i,* 2 ) 

L-j ^ o 

*1 + *2 


2*1*2 


0 


2 2 
*2 ~ x i 

-2*1*2 

0 


O' 

0 

0 - 


B 


3. On constate que la matrice de L est symetrique. Par consequent, sa decomposition 
en un tenseur asymetrique D. et un tenseur symetrique D est immediate : 

-tenseur des taux de rotation : D. — 0 

-tenseur des taux de deformation euleriens : D —L 


4. la masse volumique du milieu a l'etat initial (t = 0) est uniforme et vaut p 0 . Pour 
calculer 1'evolution de cette masse volumique, on peut evaluer le jacobien de la 
transformation en tout point et a tout instant. Dans un probleme pose en notation 
eulerienne, il est plus simple de passer par la divergence du champ de vitesse : 

dPi dV 2 dV 3 

div V — —-1- —-1- —— 

OX^ OX2 (/X3 

On sait que V 3 est constant et nul, et on a par ailleurs demontre a la question 
precedente que : 
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3V 1 _ dV 2 
dx ± dx 2 

Par consequent on a divV — 0. Un mouvement pour lequel la divergence de la vitesse 
est nulle en tout point et a tout instant est un mouvement isochore, pour lequel la 
masse volumique est constante. Par consequent, en tout point et a tout instant, on a : 

P = Po 

Le mouvement est plan, permanent, et isochore. 


5. L'acceleration y(x 1 ,x 2 ,x 3 ,t) est la derivee particulaire de la vitesse, et se calcule 
par la formule de cours suivante, qui fait apparaitre la derivee eulerienne : 

Y(x, t) = — = — + (grad v)v 

On est en mouvement permanent, done la derivee eulerienne de la vitesse est nulle : 

dV _ 

~di = ° 

Le gradient de vitesse L = grad V a deja ete calcule a la question 2, il suffit done d'en 
faire le produit contracts avec le vecteur vitesse : 


y(x, t ) = (grad V)V 


x\ + x\ 


2x ± x 2 x 2 — x\ 0" 
x 2 — x\ —2 x ± x 2 0 

0 0 0 - 


■v4(xi,x 2 ) ■ 

—x 2 - 

Xi 

B 

L 0 -1 


f(x, t ) 


A 2 (x 1 ,x 2 ) 
xl + x\ 


x 2 ■ (2x x x 2 ) + x ± ■ {x\ - xl ) 
x 2 ■ (xl - xl) + x t ■ ( 2x 1 x 2 ) 


y(x, t) 


A 2 (x v x 2 ) 

xl + xl 


-* 1*2 

-x 2 xl 

0 


x\ 




= A 2 (x 1 ,x 2 ) 


-X 1 
- X , 


6. On cherche a se placer en coordonnees polaires pour simplifier les expressions. Pour 
cela, on utilise les relations : 


(x t = r ■ cosO 
(x 2 = r ■ sinO 
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On doit simplement remplacer x t et x 2 par leurs expressions polaires dans l'expression 
de ['acceleration. On calcule d'abord la fonction A : 


A{x lt x 2 ) = 


a 


a 


a 


2n x\ + x\ 2n (r ■ cos0) 2 + (r ■ sinQ ) 2 2nr 2 


Par ailleurs on a montre dans la question precedente que l'acceleration est exprimee 
par : 


?(x, t) = - A 2 (x 1 ,x 2 ) ■ (x x ■ + x 2 ■ e^) 

y(x, t ) = —A 2 (x 1 , x 2 ) ■ (r ■ cos0 ■ + r ■ sin0 ■ e^) 

Or on a la formule de changement de repere suivante : ep = ~elcosd + ~e^sind 
Finalement: 


(X 

f(x, t ) = -A 2 0 1; x 2 ) ■ r ■ T r = - ■ ep 

L'acceleration est done purement centripete. 

7. Le champ de vitesse est donne dans l'enonce sous forme eulerienne par : 

V = -A(x lt x 2 )x 2 ■e^ + A(x 1 ,x 2 )x 1 -T 2 = A{x x ,x 2 ) ■ (-x 2 ■ e^ + ^ ■ e^) 


En coordonnees polaires, on a done : 

V = i4(;c lf x 2 ) ■ (— r ■ sin0 ■ + r ■ cos0 ■ e^) = A(jx x ,x 2 ) -r - eg 




La vitesse est purement orthoradiale (normale au rayon), et on en deduit que les 
particules ont une trajectoire circulaire autour de l'origine. Sans calcul, on obtient la 
representation lagrangienne : 


r r — R 


a 


< 

V 


e = 0 + 


2tct 2 


X 3 ~ %3 


t 


8. Pour passer des coordonnees polaires aux coordonnees cartesiennes, on utilise les 
formules suivantes : 
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, , x i X 2 

r — x\ + x? ; cos6 = — ; sind - — 

^ r r 


rz : X x X 2 

X} + Xj ; cos 0 — — ; sinQ — — 

12 R R 


On a par ailleurs Q — 


0 H--t, et done : 

2nr 2 


cosQ = cos (0 H- -t] — 

\ 27rr z ' 


et done, en developpant le cosinus : 


x 1 / cl \ / cr \ 

— - cos0 ■ cos -—- 1) — sin0 ■ sin -—- 1) 
r \/.Trr 2 / y/.Trr 2 / 


a \ a 

— = cos0 ■ cos —^- w-t — sin0 ■ sin ——^1 

r \2n(xl + x\) J \2n(xl + x\) 


— cos0 ■ cos(A(x 1 , x 2 ) ■ t) — sin0 ■ sin(y4(x 1; x 2 ) ■ t) 


x X X 

— = gt ■ cos(v4(x lf x 2 )-t)-~^- sinG4(x lf x 2 ) ■ t) 
r k k 


Finalement, on a: 


x 1 —X 1 - cos(ojt) — X 2 ■ sin(ojt) 


avec co — v4(x 1 ,x 2 ). On peut developper le sinus de la meme maniere, et obtenir 
l'expression lagrangienne de x 2 recherchee : 


/ a \ x 2 

— sin (0 H-- 1 ) = sin(0 + cot) = — 

v 2nr z ' r 


x XX 

— = sin0 ■ cos(ojt) + cos0 ■ sin(ojt) = — cos(ojt) H-sin(ojt) 

r R R 


Finalement, on a : 


x 2 — X 2 - cos(ojt) + X 1 ■ sin(ojt) 


On obtient finalement le systeme recherche : 
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' x x — X x ■ cos(ajt) — X 2 ■ sin(wt) 
■ x 2 = X 1 ■ sin(ojt) + X 2 ■ cos(ojt) 
x 3 = X 3 


9. Les trajectoires des particules sont planes et circulaires, avec une vitesse angulaire 
oo constante pour chaque particule et qui decroit en relation hyperbolique avec la 
distance a l'axe origine. L'acceleration est purement centripete. On appelle ce 
mouvement un tourbillon. 


10. On a montre a la question 3 que le tenseur des taux de deformation euleriens est 
egal au tenseur gradient de vitesse, dont la matrice dans la base B = (iq, e^, e^) est: 


= _ Afa.Xz) 

Li ^ o 

x( + 


■ 2x t x 2 
x 2 — x 1 

L o 


2 2 
X 2 — X-l 

-2x ± x 2 

0 


O' 

0 =D 

0 B 


On se place au point x — r ■ e^, et on a done x t — 0 et x 2 = r. Par consequent, la 
matrice de D exprimee dans la base B en ce point est donnee par : 


D = 


2 l( 0 , r) 


0 r 2 0 
r 2 0 0 

.0 0 0 


a 


2nr 2 


0 1 0 
10 0 
0 0 0 


On reconnait un glissement simple tel que defini dans le cours. 

La base propre d'une telle deformation est nommee B’ = (d 1; d 2 , d 3 ), avec : 


(— V2^ V2^ 

d\ - + _ 2 _e2 

< —* V2^ V2^ 
d 2 - e i - ~y 62 

v % = 

Dans cette base, la matrice de D vaut: 


D = 


a 


2nr 2 


0 0 

-1 0 

0 0 




22 



Mecanique des Milieux Continus 


Polytech Grenoble, Geo3 


3. Deformations 


Exercice A. Soit un milieu soumis a un tenseur de deformation E, dont la matrice dans 
une base B — (i e1, e 2 , e^) donnee est: 


E = 


3 

2 

A 


2 

3 

4 



1. Calculer la trace de E, ainsi que son determinant. Calculer la matrice du tenseur E 2 , 
et en deduire le deuxieme invariant du tenseur donne par l'expression classique 

E n=l{( tr E) — tr (E 2 )^. En deduire sans calcul le polynome caracteristique du 

tenseur E. 


2. Developper le determinant du tenseur E — AI et retrouver l'expression du polynome 
caracteristique. 

3. Developper differemment det[E — AT) pour l'exprimer sous la forme suivante, 
beaucoup plus interessante : det(E — AT) — —(A — E^){A — E 2 )(A — E 3 ). En deduire 

les valeurs des deformations principals et la forme de la matrice de E dans sa base 
principale. 

4. Ordonner les deformations principals, puis calculer les coordonnees du vecteur 
propre b 2 correspondant a la deformation principale intermediaire E 2 . 
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Corriges 


Exercice A. 

1. La trace de E s'obtient immediatement: 


tr E — tr 



2 

3 

4 



3 + 3 + l = 7 = £ / 


Le determinant d'une matrice 3*3 se calcule par la formule classique : somme des produits 
des termes selon les diagonales descendantes moins somme des produits des termes selon les 
diagonales montantes : 


a b c 
d e f 
g h i 


a-e-i + b- f- g+ c- d- h — g-e-c — h-f-a — i-d-b 


Applique a E, on a : 


3 2 4 

2 3 4 

4 4 1 


3*3*l + 2*4*4 + 4*2*4 — 4*3*4 — 4*3*4 — 1*2*2 


det E — —27 = E IU 

Le tenseur E 2 est le produit contracts de £ par lui-meme. Sous forme matricielle dans la 
base B — (e^e^, e^), on a : 


3 

2 

4' 


3 

2 

4' 


29 

28 

24' 

2 

3 

4 


2 

3 

4 

= 

28 

29 

24 

.4 

4 

1 . 

B 

.4 

4 

1 . 

B 

.24 

24 

33. 


La trace de ce tenseur vaut tr E 2 — 29 + 29 + 33 = 91 


Le second invariant principal de E s'obtient directement par la formule de l'enonce : 

tr !) 2 - tr (I 2 )) = i(7 2 - 91) = -21 

Par application de la formule du cours, on en deduit le polynome caracteristique du 
tenseur des deformations linearisees : 

det(E -Xl) = -A 3 + E I ■ A 2 - E u ■ A + E UI = -A 3 + 7 ■ A 2 + 21 ■ A - 27 
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2. On cherche le determinant du tenseur E — XI qui s'exprime sous forme matricielle 
dans la base B = (i e1, e^) par : 


3 

2 

4" 


1 

0 

O' 


3-A 

2 

4 ' 

2 

3 

4 

-A- 

0 

1 

0 

_ 

2 

3-A 

4 

.4 

4 

1 . 

B 

.0 

0 

1 . 

B 

. 4 

4 

1-A. 


Done : 


det^E — eT) 


3-A 2 4 

2 3-A 4 

4 4 1-A 


det(E — Al) — (3 — A) 2 *(1 — /l) + 2*4*4 + 4*2*4... 
... — 4 * (3 — /l) * 4 — 4 * 4 * (3 - 2.) — (1 — '0 * 2 * 2 


det{E -Al) = (A 2 -6A + 9)(1 - X) + 32 + 32 - 48 + 16 A -48 + 16/1-4 + 4/1 
det(E - AT) - A 2 - 6A + 9 - A 3 + 6A 2 - 92. - 36 + 36A 

det(E - Al) = -+ 3 + 7 ■ A 2 + 21 ■ A - 27 
On retrouve done par le calcul detaille le resultat de la question precedente. 

3. On cherche a obtenir un developpement du polynome caracteristique faisant intervenir 
directement ses racines, sous la forme det(E — AT) — —(2. — £’i)(2. — £” 2 ) — £ 3 )- Cette 
operation est peu aisee si on demarre de l’expression developpee du polynome, et il est 
plus judicieux de chercher cette expression des le calcul du determinant. On prend 
done comme point de depart: 


det^E — 2.7) 


3-A 2 1 

2 3-A 3 

4 4 1-A 


Selon les formules classiques de calcul de determinant, il est possible de remplacer une 
colonne (respectivement une ligne) par une combinaison lineaire des differentes 
colonnes (respectivement des differentes lignes). On remplace done la colonne Cl par 
la combinaison C1-C2 : 


3-A 2 4 


(3-A)-2 2 4 


1-A 2 4 

2 3-A 4 

— 

2 - (3 -A) 3-A 4 

— 

A-l 3-A 4 

4 4 1-A 


4-4 4 1-A 


0 

1 

>-> 
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II est egalement autorise de "sortir" une constante multiplicative d'une ligne ou d'une colonne 
du determinant, done on peut simplifier la premiere colonne : 


det(E - Al) = (A — 1) 


-1 

1 

0 


2 4 

3-A 4 
4 1-A 


On remplace ensuite la ligne LI par la combinaison lineaire L1+L2 : 


det(E - Al) = (A - 1) 


0 

1 

0 


5 -A 
3-A 
4 


8 

4 

1-A 


On a le droit de permuter deux lignes (ou bien deux colonnes) en changeant le signe du 
determinant, on peut done ecrire : 


det^E — Al) 


-a-i) 


3 - A 
3 - A 
4 


4 

8 

1-A 


Enfin, on sait qu'un determinant 3*3 de la forme obtenu peut etre simplifie en un 
determinant 2*2 par la formule suivante : 

lab 
0 c d 
0 e f 

On peut done ecrire : 

det[E -Al) = -(A - 1) 


c d 

e f 


II reste quelques calculs, mais on sait qu'on a atteint le but principal, car un 
determinant 2*2 de cette forme est un polynome de degre 2 dont on sait tres bien 
retrouver les racines. On a : 


det(E - Al) = -(A - 1)[(5 - A)( 1 -A)- 4*8] 
det{E — A) - -(A - 1)(A 2 -6 A- 27) 

Le discriminant du polynome de second degre est: 

A= (—6) 2 - 4 * 1 * (-27) = 144 
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Et les deux racines de ce polynome sont: 


6 ± VA 

2 


3 ± 


Vl44 

2 


Finalement, le polynome caracteristique du tenseur des deformations linearisees est 
donne par : 


det(E -Al) = -(A - 1 )(A - 9) (A + 3) 


Les deformations principales sont les racines de ce polynome, et on peut ecrire la 
matrice de E dans sa base principale : 



1 

0 

0 ' 

£ - 

0 

9 

0 


.0 

0 

-3. 


J Base principale 


4. On ordonne les deformations principales dans l'ordre decroissant, et on a : 


avec : 



0 

E 2 

0 



Base principale 



On pourrait calculer les coordonnees dans B — (eq, e^, e^) des trois vecteurs de la base 
principale de deformation {b ll b 2 ,b 2 \ On va le faire uniquement pour b 2 car la 
methode est toujours la meme. 


D'apres le cours, la direction principale de deformation b 2 est definie par : 


Eb 2 = E 2 b 2 

En developpant de systeme, on obtient: 


'3 2 4' 
2 3 4 


b 2 ,i 

b 2 , 2 

= E 2 

b 2 ,i 

b 2 , 2 

— 

b 2 ,i 

b 2 , 2 

.4 4 1. 

B 

P 2 ,3. 

B 

P 2 ,3. 

B 

p2,3_ 
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Dans cette expression, on a fait apparaitre les trois coordonnees (b 2: i, ^ 2 , 2 > ^ 2 , 3 ) du vecteur 

b 2 dans la base B — (e^,e^,e^). En mettant ce systeme matriciel sous forme d'un 
systeme d'equations lineaires, on a : 

3 * b 2 1 4 " 2 * b 2 2 3 ~ 4 * b 2 3 — b 2 1 
2 * b 2 ^ + 3 * b 2 2 4 " 4 * b 2 3 = b 2 2 
4 * b 2 1 3 ~ 4 * b 2 2 4 " b 2 3 = b 2 3 

Soit: 

2 * b 2 1 3 ~ 2 * /?2 2 d - 4 * /?2 3 — 0 
2 * ^21 + 2 * ^2 2 + 4 * ^2 3 — 0 
4 * /?2 ^ 3 ” 4 * /?2 2 = 0 

On remarque immediatement que b 2 ,i — ~b 22 et que b 2 3 — 0. Or on sait que le vecteur 
b 2 est unitaire (de norme egale a 1), done on peut ecrire : 

r 

b>2,\ 

< 

b>2,2 ~ 

< b 2 ,3 

On note une indetermination dans le signe de b 21 et b 2 2 . Tous les signes sont possibles, a 
condition que l'on verifie b 21 = —b 2 2 . 

Ce systeme etait trivial, mais ce n'est pas toujours vrai. Dans le cas contraire, on peut 
utiliser le pivot de Gauss, ou passer par une inversion de matrice. 
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4. Contraintes 


Exercice A. Soit un milieu continu dont l'etat de contrainte est homogene et 
stationnaire (le tenseur de Cauchy est constant dans l'espace et le temps). On considere 
au sein de ce milieu un domaine materiel D , en forme de tetraedre tel que represente 
dans la figure ci-dessous dans un repere cartesien orthonorme (0,fq, e^, e^). Les aretes 
paralleles aux axes de la base ont une longueur donnee c. 

Un experimentateur realise quelques mesures, et constate que : 

-La composante normale de la force exercee sur la face arriere (correspondant au plan 
x ± — 0) est egale a F ± 

-La composante normale de la force exercee sur la face de gauche (correspondant au 
plan x 2 — 0 ) est egale a F 2 

-La composante normale de la force exercee sur la face inferieure (correspondant au 
plan X 3 = 0) est egale a F 3 

—> _ ^ ^ ^ 

-La force exercee sur la face inclinee du tetraedre est egale a G — G 1 e 1 + G 2 e 2 + G 3 e 3 


Determiner les composantes de la matrice du tenseur de Cauchy dans la base 
( e i> e 2> e s)- 


(0,0,c) 



(0,c,0) 


Exercice B. On considere un petit cube de centre x dans un milieu continu soumis a 
des contraintes. On effectue trois experiences de chargement (a), (b), et (c), 
respectivement caracterisees par les vecteurs contraintes suivants (<r 0 etant une 
constante) : 

(a) . T(x,el) = 0 ; T(x,e^) = a 0 e^ ; T(x,e^) = a 0 e^ 

(b) . T{x,~el) = o- 0 e^ ; T(x,e^) = 0 ; T(x,e 3 ) = a Q el 


(c). T{x,el) = a 0 e^ ; f{x,e^) = a Q el ; T(x,e 3 ) = 0 
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1. Exprimer les trois tenseurs de contraintes aa (h) , et a (c) correspondant a ces trois 
cas de chargement. 

2. On effectue les trois experiences simultanement en superposant les trois systemes de 
forces. Exprimer le tenseur de contraintes correspondant a cette nouvelle experience. 
Pour la suite de l'exercice on considerera uniquement ce nouvel etat de contrainte. 

3. Calculer le vecteur contrainte exerce sur une petite surface normale a la direction n 
definie par l'expression n = eq + eq + eq. 

4. Calculer le vecteur contrainte exerce sur une petite surface de direction u normale a 
n. 


5. Donner sans calcul supplementaire les contraintes principales et les directions 
principales de contraintes. 


6. Tracer le tricercle de Mohr correspondant a cet etat de contraintes et donner la 
valeur de la contrainte de cisaillement maximale. Donner egalement la valeur de la 
contrainte moyenne et la matrice du deviateur de contraintes. 

Exercice C. En un point M d'un milieu continu et a l'instant t, la matrice du tenseur de 
Cauchy dans une base cartesienne orthonormee fixe B = (eq, iq, <qT) est donnee par 
l'expression suivante, ou a est une constante adimensionnelle : 


a = 


0.7a 

3.6a 


0 


3.6a 

2.8a 

0 


0 

0 

7.6. 


1. Montrer que le calcul des trois contraintes — 1,2,3 et l'utilisation des 

proprietes des cercles de Mohr permettent de calculer les contraintes principales qu'on 
notera cr 1 , a 2 , et <r 3 . illustrer le cas a = 1. 


2. Determiner les valeurs de a correspondant a un etat triaxial de revolution. 

3. Pour la suite de l'exercice, on pose a = 1. Determiner les directions principales de a 
en M. 


y 3 1 

4. Calculer la contrainte pour une facette de normale — —~e^ + -'e^. Tracer le point 
correspondant sur le cercle de Mohr. 

5. Determiner pour a = 1, puis pour a = 2, la valeur de la contrainte de cisaillement 
maximum ainsi que la direction de la normale correspondante. 
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Corriges 


Exercice A. 

On recherche les termes de la matrice de <r a partir d'un certain nombre de mesures de 
forces effectuees sur les faces d'un domaine materiel. L'etat de contrainte est suppose 
homogene et stationnaire, et le tenseur de Cauchy est toujours symetrique, on cherche 
done en tout et pour tout six inconnues scalaires : a ±1 , cr 22 , <r 33 , a 12 , a 13 , et <r 23 . 

On va d'abord travailler sur la force appliquee sur la face arriere du tetraedre (plan 
x ± — 0). L'enonce dit que la composante normale a cette force est egale a F 1 . On peut 
egalement remarquer que la facette en question (face triangulaire arriere) est de 
normale sortante —e[. 


On va chercher a exprimer le vecteur contrainte applique sur cette facette. On peut 
l'enoncer ainsi a partir de la formule du cours : 

f(-e^) = a ■ (—ej*) 


Avec : e ± = 


Done on a: 



eta — 


(Til 

a 12 

- a 13 


Ti-TJ = 


<T 12 °"l3 

°22 °23 


°23 

a 33. 

B 





(Til 

a 12 

013 


-1 


— 011 

<t 12 

a 22 

a 23 


0 

= 

— 012 

.<*13 

a 23 

033- 

B 

. 0 . 

B 

— 013. 


Par ailleurs, on sait que la composante normale de la force appliquee a la facette vaut 
F ± . Si on note S la surface de la facette triangulaire, on peut integrer la contrainte 
r(-e^) sur toute la surface de la facette (pour obtenir une force), et considerer sa 
composante normale (par projection sur la normale, c'est-a-dire par produit scalaire 
avec le vecteur normal sortant de la facette). On a done : 


F, = S x (f (-eD ■ (-?T)) = 5 



1 

£ 


-1 

)= S 

a 12 


0 


—013. 

B 

. 0 . 


= Sa 


11 


La surface S de la facette est donnee par S = et on peut finalement ecrire : 


2F 1 


ai1 ~ 
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Par des raisonnements identiques sur les facettes de normales —et —e^, on obtient 
de la meme maniere : 


__ 2 ^2 

a 22 — „2 

C z 

_ 2 ^*3 

<*33 ~ 2 

C z 


On a determine les trois termes diagonaux de la matrice de a. II nous reste a 
determiner les termes non-diagonaux. On va done s'interesser a la facette restante. 


Sa normale sortante est un vecteur unitaire, colineaire a la direction 


d'un vecteur d'allure generale 


rdi 

a 

LaJ 


. II s'agit done 


, ou a est un reel a determiner. Par ailleurs, ce vecteur 


doit etre unitaire, done sa norme doit etre egale a 1. On va done avoir : 


V a 2 + a 2 + a 2 — 1 


et done : 


1 



La normale sortante a la facette est done le vecteur n 


'1/V3 

1/V3 

-1/V3 


B 


On a egalement besoin de la surface S' de cette facette. II s'agit d'un triangle equilateral 
de cote cV2, et on obtient apres quelques calculs trigonometriques et/ou Pythagore : 


c 2 V3 


On peut maintenant ecrire le vecteur contrainte subi par cette facette en n'importe 
lequel de ses points (puisque l'etat de contrainte est homogene): 

T(n ) - an 


et done : 
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?(n) = 


(Til 

a Y2 

a 13 


'1/V3' 

a 12 

°22 

a 23 


1/V3 

Pl3 

a 23 

a 33. 

B 

-1/V3. 


T(n) = 


V3 


(Til + a l2 + °13 
cr 12 + °22 + °23 
013 + ff 23 + ff 33 


L'enonce precise que la force subie par la facette est mesuree comme etant egale a : 


G — G 1 e 1 + G 2 e 2 + G 3 e 3 

Or cette force est egalement le produit du vecteur contrainte (homogene) par la surface 
de la facette, on a done : 


°11 + *-*12 + < Ti3 

G — G^e^ + G 2 6 2 + G 3 e 3 — S r • T(n ) = —= <Ti 2 + <r 3 2 + 023 

^ -013 + ff 23 + <T 33 

Cette equation matricielle nous fournit trois equations scalaires : 


( S' 

G t = — (a lt + <r 12 + <j 13 ) 
S' 

~ (°"l2 + cr 22 + <T 23 ) 

5' 

^3 = (°"l3 + °"23 + ^ 33 ) 


On peut alors remplacer les termes diagonaux de la matrice par leurs valeurs trouvees 
en debut d'exercice, et faire de meme avec la surface de la facette : 


r 

Gi 



c 2 (2F 1 \ 

c 2 ( 2F 2 \ 

2 ( ff 12 + C 2 + °23 j 

c 2 ( 2 F 3 \ 

2 r 13 + ff23 + c 2 ) 


On dispose alors d'un systeme de trois equations a trois inconnues : <Ti 2 , a 13 , et <r 23 . 
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r 2 G 1 — 2 F 1 

ff 12 + ^13 = ~ 2 

2 G 2 -2 F 2 

a 12 + a 23 — 

c z 

2G 3 -2F 3 

,°"l3 + a 23 - ~z 

On obtient finalement les resultats demandes suivants : 



(gi ~ Fi) + (gz ~ f 2 ) - (£3 ~ f 3 ) 

c 2 

(gi ~ FQ - (G 2 - F 2 ) + (G 3 - F 3 ) 
c 2 

-(G x - FO + (G 2 - F 2 ) + (C 3 - F 3 ) 


Exercice B. 

1. Considerons l'etat de chargement (a). On nous donne trois equations vectorielles : 
T(x,el) = 0 ; T(x,e^) = a Q T 3 ; T(x,e^) = a Q T 2 

D'apres la premiere equation, on peut ecrire : 

"0n <t 12 cr 13 “i ru r^n" 

T(x, e^) - ■ eq = 012 <r 22 <r 23 ■ 0 — 012 = 0 

013 023 033 J LoJ L°13- 

On en deduit directement que — cr 12 — a 13 — 0 

On considere maintenant la deuxieme equation, qui nous permet d'ecrire : 

|"0n <t 12 a 13 i mi [& 12 

T(x,e£) = (7 (a) ■ el = 012 022 023 ■ 1 = 022 = a Q ~el 

.013 023 033J LOJ L023. 

On en deduit que <r 12 = a 22 — 0 mais aussi que <r 23 = a 0 . 

Enfin, on considere la troisieme equation, qui donne directement: 

■0-li cr 12 a 13 0" "0i3 _ 

T{x, tQ = ■ ~el = 012 022 023 ■ 0 = 023 - cr 0 e^ 

-013 023 033-1 LlJ 1-033- 

On en deduit que <r 13 = a 33 = 0 et on confirme que o 23 — a 0 . 
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On a done calcule les six termes independants de la matrice du tenseur de Cauchy dans 
la base B — (eq, e 2 , e^), et on peut ecrire cette matrice : 



0 

0 

.0 


0 

0 

0 O 


0 - 

0 O 


On rappelle que cette matrice est uniquement valable pour l'etat de chargement note 
(a) dans l'enonce. En appliquant des raisonnements identiques, on obtient egalement 
les matrices des tenseurs de contraintes correspondant aux etats de chargement (b) et 
(c) : 



0 

0 

-°o 


0 

0 

0 


0 O 

0 


a (c) - 


0 

0 O 

0 


0 O 

0 

0 


2. Sans calcul supplementaire, on peut tirer parti du fait que la superposition de 
plusieurs etats de contraintes correspond a la somme des tenseurs de contraintes 
correspondants. Si on superpose les chargements (a), (b), et (c), on obtient done le 
tenseur de contraintes dont la matrice dans la base B = (e^e^, e^) s'exprime par : 


0(a)+(O + (c) - 0(a) + 0(b) + 0(c) - 


0 

0o 

On 


00 

0 


00 

<Jn 


<T, 


0 


3. On cherche a exprimer le vecteur contrainte sur une facette dont on nous donne la 
direction n. 11 nous faut d'abord determiner les coordonnees du vecteur unitaire normal 
a cette facette. Ce vecteur est colineaire a n, mais de norme egale a 1. On peut le noter 

n' et l'exprimer ainsi : 


—j n n 1 

n = 

Le vecteur contrainte applique a la facette definie par cette normale unitaire est ensuite 


donne directement par la formule du cours : 




_ 0 

0o 

0o" 

1 [1 

T(x, n) = T ( x, n') = a {a)+(b)+(c) ■ n' = a 0 

-0o 

0 

0o 

0o 

'vf 1 

L1J b 


Et done : 


T 

1 
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T ( x, n ) 


2a 0 

V3 



2a 0 

V3 


ft = 2<r 0 n' 


4. On cherche a determiner le vecteur contrainte applique a une facette de normale it, 
dont on sait seulement que cette normale est egalement normale au vecteur n defini 
plus haut. On cherche done T(x, it), sachant que it ■ n = 0. Ecrivons ces deux donnees 
sous forme mathematique : 


it ■ it = + u 2 n 2 + u 3 n 3 = u ± + u 2 + u 3 — 0 


T ( x, it) 


a (a)+(b)+(c) 



0 

00 

0o" 


~U ± ' 


U 2 + U 3 

it = 

0o 

0 

00 


u 2 

= 00 

IL ± + U 3 


-0o 

00 

0. 

B 

U 3 _ 

B 

U ± + U 2 _ 


On en deduit: 


T(x,u ) = — a 0 


u t 

u 2 

L u 3i 


— —a 0 u 


5. D'apres les formules precedentes, on peut formuler plusieurs observations. On a 
montre que : 

T (jc, n'j = 2<j 0 n' 

Par consequent la direction n' est une des directions principales de contrainte. Le 

vecteur unitaire n' est done l'un des vecteurs de la base principale de contraintes, et la 
contrainte principale correspondante est egale a 2<r 0 . 

Par ailleurs, on a montre que, pour tout vecteur it perpendiculaire a n ', on a : 


T(x,u ) = — a 0 u 


On en deduit que —a 0 est une racine double et que le tenseur de Cauchy est 
cylindrique. Dans sa base principale, sa matrice s'exprime par : 


<7i 


(a)+0)+(c) 


2ao 
0 
L 0 


0 

-0o 

0 


0 - 
0 


-OoJ 


Base principale 


La base principale est formee du vecteur n' et de tout couple de vecteurs 
perpendiculaires appartenant au plan normal a n'. 
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6. En ordonnant les valeurs de contraintes (et en supposant que <r 0 > 0), on a : 


<T (a)+(fo)+(c) 


<7i 

0 

0- 

0 

^2 

0 

O 

0 

^3- 


Base principale 


avec : a ± — 2a 0 et a 2 — °3 — ~ a o- Le tricercle de Mohr a Failure suivante : 



Les valeurs propres <r 2 et a 3 sont confondues, done le tricercle se resume a un cercle 
unique. Par ailleurs, la contrainte de cisaillement maximale et la contrainte moyenne 
sont donnees par les formules de cours : 


T-max 



1 _ <Ti + a 2 + (To 

0Vn = 3 tra = ---= 0 


La contrainte moyenne est nulle, done le tenseur de contraintes n'a pas de partie 
spherique, et le deviateur de a est egal a a. 


Exercice C. 

1. On dispose de la matrice du tenseur des contraintes dans la base B — (eq, e^, (%) : 


a — 


"0.7 a 
3.6a 


L 0 


3.6a 

2.8a 

0 


0 " 
0 

7.6. 
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Suivons l'enonce, et calculons les trois vecteurs contraintes existant au point M, pour 
des facettes de normales respectives el, el, et ~el : 



0.7a 3.6a 0 " 


1 


"0.7a" 

T(el) = ael = 

3.6a 2.8a 0 


0 

= 

3.6a 


. 0 0 7.6. 


.0. 


. 0 . 


"0.7a 3.6a 0 


"0 


"3.6a" 

T(el) = ael = 

3.6a 2.8a 0 


1 


2.8a 


. 0 0 7.6. 


.0. 


. 0 . 


0.7 del + 3.6 del 


3.6 del + 2.8ae^ 


T(el) = a el 


"0.7a 3.6a 0 " 


'O' 


■ 0 " 

3.6a 2.8a 0 


0 

= 

0 

o 

o 

CTN 


.1. 


.7.6. 


7.6el 


Chacun de ces vecteurs peut etre decompose en une composante normale et une 
composante tangentielle, afin d'etre representes dans le plan de Mohr. Par exemple, 
pour le vecteur contrainte applique a la facette de normale el, on peut ecrire : 


T(el) = T nl -el + T tl -tl 


Dans cette equation tres generale, t x est un vecteur unitaire appartenant au plan de la 
facette, et T n i et T n sont respectivement la composante normale et la composante 
tangentielle du vecteur contrainte. Par identification avec l'expression trouvee plus 
haut, on peut ecrire : 


h 



T n i - 0.7a et T tl — 3.6 a 

En appli quant la meme operation a T(el) et T(el ), on peut ecrire de la me me maniere 


T n2 = 2.8a et T t2 = 3.6a 
T n 3 — 7.6 et T t3 — 0 


Avec ces resultats, on peut tout a fait placer sur le plan de Mohr les points 
correspondant a chacun de ces trois vecteurs contraintes. 11 s'agira des points de 
coordonnees (T nl ,T tl ), (T n2 ,T t2 ), et ( T n3 ,T t3 ). 11 faut noter pour l'instant que la 
connaissance de ces points n'est pas suffisante pour tracer le tricercle de Mohr et done 
pour determiner les contraintes principales. 

On doit maintenant appliquer un certain nombre de raisonnements lies a des proprietes 
du cercle de Mohr. On remarque que la composante tangentielle de T (el) est nulle, ce 
qui signifie par definition que el est une direction principal de contrainte. On en deduit 
done que les vecteurs el et el sont situes dans le plan des deux autres directions 
principales (pour l'instant inconnues), et que les points du plan de Mohr correspondant 


38 



Mecanique des Milieux Continus 


Polytech Grenoble, Geo3 


aux deux vecteurs contraintes lies a ces facettes sont situees sur le meme cercle de 
Mohr (c'est la definition d'un cercle donne dans le cours). Par consequent, les points 
(T nl , T tl ) et (T n2 ,T t2 ) appartiennent tous les deux au meme cercle de Mohr. C'est 
aussi le cas de leurs symetriques par rapport a l'axe horizontal, c'est-a-dire des points 
(T n i, —T t i) et (T n2 , —T t2 ). 


On tire enfin parti des regies d'orientation donnees dans le cours, et on observe que les 
deux facettes de normales ~el et sont orthogonales, et done que les points qui leur 
correspondent dans le plan de Mohr sont necessairement diametralement opposes sur 
le cercle de Mohr auquel ils appartiennent. On en deduit que les coordonnees du centre 
du cercle de Mohr que l'on recherche sont au milieu d'un segment reliant deux points 
diametralement opposes : 



+ T n 2 ^ (0.7a + 2.8a 

2 ’°) ~ V 2 



(1.75a, 0) 


Le rayon de ce cercle se calcule ensuite par Pythagore : 

R = 7(2.8a - 1.75a) 2 + (3.6a) 2 = aV(2.8 - 1.75) 2 + (3.6) 2 = 3.75a 


On en deduit qu'une des contraintes principales est egale a 1.75a — 3.75a = —2a, et 
qu'une autre contrainte principale vaut 1.75a + 3.75a = 5.5a. Quant a la troisieme, 
on l'a deja determinee implicitement au debut de l'exercice, et elle vaut 7.6. La 
situation est representee dans le plan de Mohr pour a = 1: 
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2. Un etat triaxial de revolution correspond a n'importe quel etat de contraintes pour 
lequel deux des contraintes principales sont egales. 

Dans notre cas, sans accorder d'attention a l'ordre des contraintes principales (puisque 
celui-ci depend de a ), on peut enoncer que les contraintes principales sont egales a : 

a ± = —2a ; o 2 — 5.5a; <r 3 = 7.6 

On peut relever dans le cas present trois situations pour lesquelles l'etat de contraintes 
est triaxial de revolution : 


<Ti = o 2 => —2a — 5.5 a => a — 0 
o"i = a 3 => —2a — 7.6 => a = —3.8 
° 2 — °3 =* 5.5a = 7.6 => a = 1.38 


3. On pose a — 1, on ordonne les contraintes, et on a done o 1 = 7.6, a 2 — 5.5, et 
a 3 = -2. 

On cherche les directions principales de contraintes (cq, c^, c^). Plus precisement, on 
cherche a trouver les coordonnees de ces trois vecteur propres dans la base d'origine 
(e 1; e 2 , e 3 ). 


On a deja demontre que ~el est une direction propre, et lorsque a — 1 cette direction 
propre est celle de la contrainte principale maximale. On peut done enoncer sans calcul 
que : 


c 1 = e, = 


On cherche alors les coordonnees des deux autres vecteurs principaux, tels que 


c 2 — c 2 ,l ' e l + c 2 2 ■ 6 2 + C 2 3 ■ e 3 


c 3 — c 3 ,i ' e i + c 3 2 ■ e 2 + c 3 3 ■ e 3 


Puisque ~cl — e^, les vecteurs et sont orthogonaux a e^, et on a c 2 3 = c 3 3 — 0. 
II nous reste quatre scalaires a trouver : c 21 , c 22 , c 31 ,etc 32 . 


On applique la formule classique definissant une direction principale de contrainte 


GCi t-i 
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0.7 3.6 0 ] rc 2;1 -| r0.7c 2(1 + 3.6c 2 2 

or on a = 3.6 2.8 0 c 2;2 = 3.6c 2 ± + 2.8c 2 , 2 

- 0 0 7.6J b L 0 J B l 0 

Done on obtient ac 2 — (0.7 c 2tl + 3.6c 2 2 )e^ + (3.6c 2 , ± + 2.8 c 22 )el 

Par ailleurs on a (Tel — a 2 ~cl = 5.5 ■ [c 2/1 e^ + c 2 2 el] 

On en deduit le systeme suivant: 

f5.5c 21 = 0.7c 2 x + 3.6 c 2 2 
(5.5c 2 2 = 3.6c 2 x + 2.8c 2 2 

r—4.8 c 2 i + 3.6c 2 2 = 0 
l 3.6 c 2( i — 2.7c 22 — 0 

-4.8 3.6 r c 2 ,ii = roi 

3.6 -2.7J \- c 2,2\ LoJ 

La solution de ce systeme est triviale : c 21 — 0 et c 2 2 — 0. Cette solution n'est 
evidemment pas celle qui nous interesse. En se penchant un peu plus sur le systeme, 
on observe que ses deux equations sont en fait identiques. C'est le genre de piege que 
les calculs peuvent faire apparaitre. 

11 y a une information que l'on n'a pas encore utilisee : le fait que le vecteur doit etre 
unitaire. On a done un nouveau systeme : 

r—4.8 c 2) i + 3.6c 2 2 — 0 
l cl ,! + C 2 2 = 1 

3 

De la premiere equation il vient c 21 = -c 2;2 , que Ton peut integrer dans la seconde 
equation : 

, 9 , 

° 2 ’ 2 + 16 C2 ' 2 ~ 1 

Finalement, on obtient: c 2 2 — 0.8 et c 21 = 0.6. On note enfin : 

0.6' 
cl = 0.8 

- 0 J B 

On pourrait appliquer la meme methode pour obtenir c 31 , et c 3 2 , mais on peut aussi 
tirer parti du fait que (c^c^c^) est une base orthonormee directe. On peut done 
utiliser le produit vectoriel et ecrire : 



41 



Mecanique des Milieux Continus 


Polytech Grenoble, Geo3 


Done on a : 


C 3 — Ci A C 2 


'O' 


'0.6' 


-0.8 

0 

A 

0.8 


0.6 

.1. 

B 

. 0 . 

B 

. 0 . 


\ 3 1 

4. On definit la facette de normale + Le vecteur contrainte subi par 

cette facette s'obtient directement par la formule du cours : 


T(n±) — 


0.7 

3.6 

0 ' 


V3/2 


'2.4' 

3.6 

2.8 

0 


1/2 

= 

4.5 

. 0 

0 

7.6. 

B 

- 0 - 

B 

. 0 . 


Pour tracer le point correspondant a ce vecteur contrainte dans le plan de Mohr, on a 
besoin d'une composante normale et d'une composante tangentielle. La composante 
normale s'obtient, comme d'habitude, par projection du vecteur T(n^) sur la direction 
normale a la facette, e'est-a-dire par un produit scalaire : 



'2.4' 


V3/2 

T n = T (n/) ■ n/ = 

4.5 


1/2 


. 0 . 

B 

- 0 - 


Le vecteur contrainte L(n/) est la somme d'un vecteur T n n/ oriente dans la direction 
normale et d'un autre vecteur oriente dans une direction tangentielle t ± qui n'est pas 
connue. La composante tangentielle est done la norme du vecteur obtenu en 
retranchant la composante normale du vecteur T(n/) a lui-meme : 


T t = || T(jh) - T n nl\ 


[2.4' 


V3/2 


4.5 

- 4.3 

1/2 


i 0 - 

B 

- 0 - 

B 


On peut alors reporter le point correspondant sur le graphique (voir figure precedente). 
On constate qu'il appartient au cercle 1. C'est logique, car ce point est defini pour une 

facette de direction normale n/ = — ej* + -e^ qui appartient au plan des deux 
directions principales (c^, c^) (c'est le plan normal a e^). 
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5. Pour a — 1, on a o x — 7.6 et a 3 — —2 . La contrainte de cisaillement maximale est 
done donnee par : 


Qi <J 3 

t = —-- = 48 

L max 2 

L'abscisse du point correspondant sur le plan de Mohr correspond a celle du centre du 
grand cercle de Mohr, soit: 


Le point du cercle de Mohr correspondant a la facette de cisaillement maximal est 
done de coordonnees (2.8,4.8) (voir figure). 

Par definition, la facette qui recoit cette contrainte possede une normale qui appartient 
au plan des deux directions principales q et (qu'on appelle "plan de cisaillement 
maximal"). En fait, la facette de cisaillement maximal correspond toujours a une 
orientation situe selon la bissectrice des directions principales q et cq, qui est done 
donnee par : 



00 

o 

1 

o 


cl + cl _ 

2 

0 + 0.6 


-0.4' 



— 

0.3 

2 

2 

1 + 0 


. 0.5 . 


9 



L 2 J 


11 faut noter que ce vecteur n'est pas unitaire (meme s'il suffit a definir la direction de 
la normale a la facette). 

Pour a = 2, la methode est strictement la meme, excepte que le vecteur cq a change, et 
il faut done le recalculer pour obtenir la bonne orientation de la facette de cisaillement 
maximal. 
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5. Elasticity 


Exercice A. Soit le massif rectangulaire de grande longueur represente sur la figure 
suivante : 



On pose plusieurs hypotheses : 

-Comportement elastique lineaire, de parametres Eel v. 

-Face superieure soumise a un pression uniforme p. 

-Face inferieure en appui glissant et non-frottant, fixee seulement sur l'axe (0,e^). 
-Forces de volume negligeables. 

-Hypothese des petites perturbations. 

-On se place dans la zone centrale du massif pour se premunir des effets de bord. 
-Champ de deplacement donne par U — (A ■ X 1 + C ■ A 2 )e^ + (—C ■ X ± + B ■ X 2 )E^. 

1. Representer la section du massif contenant le point M, ainsi que les conditions aux 
limites sur cette section. 

2. Calculer le tenseur des deformations s, et donner ses valeurs principals et ses 
directions principals. De quel type de deformation s'agit-il ? 

3. Exprimer analytiquement le tenseur de Cauchy a (utiliser les coefficients de Lame 
pour simplifier les notations). 

4. En utilisant les conditions aux limites en contraintes, calculer les expressions des 
parametres A et B. En utilisant les conditions aux limites en deplacement, calculer C. 

5. On pose E — 210000 MPa et v = 0.3 (Massif en acier). Les dimensions de la 
section sont L — 0.2m et H — 0.15 m, et le chargement vaut p = 500 MPa. Calculer 
le vecteur deplacement du point M. Tracer la deformee de la section. 
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Exercice B. On s'interesse a un barrage poids de section triangulaire, dont la geometrie 
est proposee sur la figure suivante. Pour rendre le probleme plan, on s'interesse en 
particulier a une "tranche" de ce barrage situee suffisamment loin des massifs 
d'ancrage lateraux. Le barrage est soumis a son poids p b g et a la poussee p e gx 2 de 
l'eau, la surface de celle-ci etant situee a une hauteur h egale a celle du barrage. A sa 
base, le barrage est parfaitement encastre. 

Le materiau composant le barrage est suppose elastique lineaire et isotrope, de 
parametres Let v. On neglige la pression atmospherique. 



1. Expliquer pourquoi il s'agit d'un probleme de deformations planes, et donner la 
forme generale des inconnues U, s, et a. 

2. Decrire precisement les conditions aux limites d'une "tranche" du barrage. Ecrire les 
equations d'equilibre mecanique, d'abord a l'interieur du domaine puis au niveau de ses 
conditions limites en contraintes. 

3. On postule que tous les termes du tenseur de contrainte sont des fonctions lineaires 
des coordonnees x ± et x 2 . Donner l'expression generale des termes de la matrice d'un 
tel champ de contrainte. Montrer qu'il verifie l'equation de Beltrami, et est done la 
solution du probleme elastique. 
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Rappel: l'equation de Beltrami est donnee par : 

1 —-♦ = = 

Act + —— grad grad tr a = 0 

Le laplacien d'un tenseur (par exemple Act) est le tenseur dont chaque terme est le 
laplacien des termes du tenseur d'origine (ici a). On a done (Act)^- = A(ct £; ). Les 
autres termes de l'equation sont calculables. 


4. A l'aide de la loi de Hooke, montrer que ct 33 = v(a 1± + ct 22 ). 

5. A l'aide des equations d'equilibre et des conditions aux limites en contraintes, 
calculer ct en tout point du systeme. A l'aide de la loi de Hooke, en deduire £ en tout 
point du systeme. 

6. Montrer que la contrainte tangentielle la plus intense (en valeur absolue) du systeme 
est obtenue au point B, et calculer cette contrainte. 
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Corriges 


Exercice A. 

1. La section du massif contenant le point M est un rectangle. II est soumis aux 
conditions aux limites suivantes : 

-face superieure : pression uniforme egale a p. 

-faces laterales : pression uniforme egale a 0. 

-face inferieure : deplacement vertical bloque, deplacement horizontal bloque 
au point central. 



2. Le champ de deplacement est donne par l'expression : 


f A ■X 1 + C-X 2 
U = -C ■X 1 + B -X 2 

i 0 


II s'agit manifestement d'un mouvement plan. Le tenseur des deformations linearisees 
se calcule par : 


£ = - (jqrad U + (grad f/) j 


Or grad U — 


1 ^ 

dU 1 

dVil 




dx 1 

dx 2 

dx 3 


' A 

c 

O' 

dU 2 

dU 2 

dU 2 




— 

— 

— 

— 

-C 

B 

0 

dx 1 

dx 2 

dx 3 



dU 3 

dU 3 

dU 3 


. 0 

0 

0 . 

-dx t 

dx 2 

dx 3 - 
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Finalement, £ — 


A 

0 

.0 


0 

B 

0 


O' 

0 

0. 


On constate qu'il s'agit d'une deformation plane, que la base principale de deformation 
est la base d'origine (ey, e^, e^), et que les deformations principales sont {A, B, 0}. 

3. On applique la loi de Hooke formulee en fonction des coefficients de Lame : 

g — 2ji ■ E + A ■ tr(J) ■ I 


avec ii — 


et A = 


Ev 


2(l+v) (l+v)(l—2v) 

On a tr(E) — A + B, et il vient done : 


Et finalement: 

a = 



A 

0 

O' 


1 

0 

O' 

a — 2/i- 

0 

B 

0 

+ A ■ (A + B) ■ 

0 

1 

0 


.0 

0 

0. 


.0 

0 

1. 


2 ft A + A • (A + B) 0 0 

0 2/iB + A - (A + B) 0 

0 0 A ■ (A + B) 


Les tenseurs a et £ sont homogenes, e'est-a-dire qu'ils sont egaux en tout point du 
sy steme. 

4. Sur la face superieure, la pression appliquee peut se noter : 


T d = ~P■ e 1 


Par ailleurs, la condition aux limites d'equilibre sur cette face superieure s'ecrit: 

on = T d 

La normale unitaire sortante a cette surface est —e[, done on a <r(—e^) = T d 
Par consequent, on peut ecrire (2 fiA + A ■ (A + 5)) ■ ~el — — p ■ ~e{ 

Et finalement: 2/iA + A ■ (A + B) = p 

On applique la meme methode sur une des faces laterales, et on obtient: 

2 fiB + A ■ (v4 + 5) = 0 


11 en resulte un systeme lineaire de deux equations a deux inconnues A et B : 
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((2p + X)A + AB> — p 
UA + (2p + A)B = 0 


La deuxieme equation donne B = 
obtenir A + A — ) = p. 


-XA 

2 jU+A 


que Ton reinjecte dans la premiere pour 


Apres calcul, il vient: 


— 2p — A A 

A=P 4 p(p + A) 6t B = V 4 p(p + A) 

Sur la face inferieure, le deplacement est bloque verticalement. En posant X 1 — 0, on 
peut done ecrire : 

U ± (0,X 2 ,X 3 ) — 0 


Or U x (X ± ,X 2 ,X 3 ) = A ■ X x + C ■ X 2 . 

On en deduit que C ■ X 2 — 0 V(A 2 , X 3 \ et done que C = 0. 


5. En utilisant les resultats precedents, on peut reecrire le champ de deplacement 


f —2u — A 

A ' Xl = P 4p(p + A)' Xl 

U = < A 

B-X 2 =p — 7 -— ■ A 2 

F 4p(p + X) 

0 


Les coordonnees du point M sont X 1 — H — 0.15m et X 2 — L — 0.2 m. 


Par ailleurs, les coefficients de Lame valent: 


210000 


B = 


A = 


2(1+ v) 2(1+ 0.3) 
210000 * 0.3 


= 80770 MPa 


= 121150 MPa 


(1 + 0.3)(1 - 2 * 0.3) 

11 vient finalement A — —2.1666 ■ 10~ 3 et B — 9.2855 ■ 10~ 4 , et enfin 


(— 0.325mm 
U (M) = | 0.186mm 

l 0 
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Exercice B. 

1. II s'agit d'un probleme de deformations planes car la section du barrage est identique 
sur toute sa longueur, de meme que le chargement qu'il subit. On en deduit que son 
champ de deplacement est entierement situe dans le plan perpendiculaire a sa longueur 
(plan d'une "tranche"), et que sa deformation dans la direction de sa longueur est 
empechee par les tranches voisines. Les champs inconnus ont la forme suivante : 



[U 1 (X 1 ,X 2 ) 

£ 11 £ 12 0 


°ii °’i2 0 

u = < 

u 2 (.x lt x 2 ) ; l = 

£ 12 £ 22 0 

; a = 

®\2 g 22 0 


0 

- 0 0 0 - 


- 0 0 <^ 33 - 


Par ailleurs, en elasticity lineaire on applique systematiquement l'HPP, et on peut 

—> 

utiliser indifferemment les vecteurs position initial X et final x. 

2. Sur une "tranche" du barrage, les conditions aux limites sont les suivantes : 

-Sur sa base il est parfaitement encastre : U — 0 sur S uAB 

-Sur la face immergee (notee S WA ), il subit la pression de l'eau 

-Sur la face a fair libre (notee S t0B ), il subit la pression atmospherique qui est 

jugee negligeable. 

Les equations d'equilibre du milieu continu sont: 

-A l'interieur du domaine, forme locale du PFD appliquee a l'equilibre statique 
du solide charge : 

div a + pg = 0 sur D 

-Sur la surface immergee, la normale sortante est n = — eq, et la contrainte 
subie par la surface est proportionnelle a la profondeur x 2 et orientee dans la 
direction de ~e[ : 

o (~eD = P e g *2 ' eT sur S t0A 

^ ^ _ > _ > 

-Sur la surface a fair libre, la normale sortante est n = (e ± — e 2 ), et la 

contrainte subie par la surface est nulle : 

= 0 sur S t0B 

3. L'enonce postule que tous les termes du tenseur de contrainte sont des fonctions 
lineaires des coordonnees x ± et x 2 , on a done proportionnalite entre chaque terme de 
ce tenseur et ces coordonnees. On peut done ecrire les relations suivantes : 

-f b 11 x 2 
*^12 ~ a 12 X l "f ^12 X 2 

&22 ~ ^22-^1 "f ^22 X 2 

@33 ~ a 33 X l + b 3 3 x 2 
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Dans ces expressions, les termes {a llt b lt , a 12 , b 12 , a 22 , b 22 , a 33 , b 33 } sont des 
constantes de proportionnalite qui sont pour l'instant inconnues. 

Quant aux deux autres composantes du tenseur de Cauchy, on a deja demontre qu'elles 
sont nulles car on est en deformations planes : a 13 — o 23 — 0. 

Pour que ce champ de contrainte soit admissible, il doit verifier l'equation de Beltrami: 


1 —-. = = 

A<t + —j -—grad grad tr a — 0 

On ecrit d'abord les differents termes de cette equation sous forme indicielle. Le 
premier terme permet d'ecrire : 


, d 2 a ii d 2 a ii d 2 a ii 
(Aff)y = A (a 1; ) = + -Gt + 1 ‘ 


dx( 


dx\ 


dx\ 


le deuxieme terme est legerement plus long a deriver. On developpe d'abord la trace 
(grad grad tr a) — (grad grad{a xl + <r 22 + a 33 )) 


On developpe ensuite le gradient du champ scalaire de cette trace, qui est done un 
champ vectoriel: 


/ 


(grad grad tr a) 


grad 




dx ± 

d 


dx 2 

d 


L dx. 


(<Tn + < 7 22 + <733) 
((Tn + <T 2 2 + O33) 


(<Tn + <J 22 + <T 33 ) 


l\ 


I/.. 

IJ 


On calcule enfin le gradient de ce champ vectoriel, qui est done un champ tensoriel 


(grad grad tr a).. 


d 2 

d 2 

d 2 

2 ^ a H + a 22 + °3 3 ) 

d Xl dx 2 (ffll + ff22+CT 33 ) 

d Xl dx 3 ( ff ll + ff 22 + T 3 3 ) 

d 2 

d 2 

d 2 

d Xl dx 2 (ail+Cr22 + CT 3s) 

Oil + °22 + ^33) 
dx 2 

dx 2 dx 3 (ffl1 + 022 + ° 33) 

d 2 

d 2 

d 2 

d Xl dx 3 Kl + C722 + C7 33 ) 

dx 2 dx 3 ^ 11 + a22 + a33 ' ) 

^2 Oil + °22 + °3 3 ) 


i. . 

1; 
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II vient: 


(,grad grad tr a) = 


dx t dXj 


(°11 + <7 22 + 033 ) — 


d 2 a 


11 


+ 


d 2 a- 


22 


+ 


d 2 a- 


33 


dx t dXj dx t dxj dx t dXj 


Finalement, l'equation de Beltrami peut s'ecrire en notation indicielle : 

' d 2 <Tn d 2 a ->9 d 2 a- 


d 2 a ii d 2 a ii d 2 a tj 1 

+ -r4- + -r4- + 


dx( ' dx 2 p ' v ' 2 


+ 


+ 


33 


dx 3 1 + v \dx t dXj dx t dXj dx t dxj 


= 0 


II faut signaler que la notation d’Einstein aurait permis de reduire sensiblement la 
longueur de cette demonstration. 


On voit que cette derniere equation ne fait apparaitre que des derivees secondes des 
termes du tenseur de contrainte par rapport aux coordonnees spatiales. Or ces termes 
suivent une variation affine par rapport aux coordonnees spatiales, done toutes ces 
derivees secondes sont done nulles. II en resulte que l’equation de Beltrami est verifiee, 
et que l’hypothese d’un champ de contrainte affine est satisfaisante. 


4. On cherche a utiliser la loi de Hooke pour demontrer que <t 33 = + a 22 )- On 

l’utilise d’abord dans le sens contraintes-deformations : 


-011- 



-1 - V 

V 

V 

0 

0 

0 - 


f ii" 

a 22 



V 

1 -v 

V 

0 

0 

0 


f 22 

CT 33 

E 


V 

V 

1 -v 

0 

0 

0 


f 33 

CT 12 

“(l + v)(l 

— 2v) 

0 

0 

0 

1 - 2v 

0 

0 


f 12 

CT 13 



0 

0 

0 

0 

1 - 2v 

0 


f 13 

-023- 



- 0 

0 

0 

0 

0 

1 - 2v- 


- f 23- 


II en resulte que fr 33 = (1+v)(1 _ 2v ; (v£n + ve 22 + (1 - v)£ 33 ) 


On sait que s 33 = 0 en deformations planes. On utilise alors la loi de Hooke dans le 
sens deformations-contraintes pour exprimer £ ±1 et £ 22 : 


£ 11" 


- 1 

—v 

—v 

0 

0 

0 - 


~ a ii~ 

£ 22 


—v 

1 

—v 

0 

0 

0 


g 22 

S 33 

1 

—v 

—v 

1 

0 

0 

0 


°33 

S 12 

“ E 

0 

0 

0 

1 + V 

0 

0 


°12 

£ 13 


0 

0 

0 

0 

1 + V 

0 


°13 

£ 23- 


- 0 

0 

0 

0 

0 

1 + V- 


- <J 23- 


1 1 

On en deduit: £ ± i = - (ffu - Vfr 22 - Vfr 33 ) et £ 22 = - (-va 11 +a 22 - vfr 33 ) 
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On reinjecte ces deux expressions dans la formule de a 33 : 


a 33 ~ 


(l + v)(l-2v)G?( ff11 “ Vff22 " VCT33) + ! ( - W " +ff22 " Vff33) 


V 


a 33 “ 


(1 + v)(l - 2v) 


(cTnCl - v) + o 22 ( 1 - v) + <t 33 (-2v)) 


On isole ensuite <r 33 , et on obtient: 

( 2v 2 

<t 33 I 1 + 


v(l-v) 


(l + v)(l-2v)/ (l+v)(l-2v) 


(°ii + ^22) 


Et finalement: 


v(l-v) 


°33 — 


(1 + v)(l - 2v) 


1 + 


2v 2 


On + a 22 ) = v(o 1± + c t 22 ) 


(1 + v)(l - 2v) 


La derniere simplification est laissee au lecteur, mais elle a ete verifiee ! 


5. Les equations d'equilibre ont deja ete donnees a la question 2, on va done 
maintenant les utiliser en tirant parti de l'expression generale des termes du tenseur de 
contraintes donnes a la question 3. L'objectif est de trouver les expressions des 
constantes {a±i, b 1± , a 12 , b 12 , a 22 , b 22 , a 33 , b 33 }, et done de definir entierement le 
tenseur de Cauchy dans le systeme. 


La premiere equation d'equilibre est celle liee a la forme locale du PLD, valable pour 
tout point du systeme : 

d iv o + pg — 0 sur D 


On calcule d'abord la divergence du tenseur de Cauchy, donnee par : 


(div <r). = 


+ do i2 do i3 

dx 1 dx 2 dx 3 


Les vecteurs div o et g s'expriment done par : 


div o = 


d CTi 

~dx^ + 
d a i 2 
dx ± 
d a i3 
_ dx ± 


do 


12 


dXn 


do- 


22 


dx 7 


do- 


23 


dxn 


+ 


+ 


+ 


do 13 



dx 3 

^^23 

et g = 

'O' 

9 

- 0 - 

dx 3 
d @33 

dx 3 \ 
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Puisque a 13 — a 23 — 0, on peut done ecrire le PFD sous la forme 


( d Oi i d(Ji 

+ 


dxi dx 7 


= 0 


d<J12 d<7 2 2 

— + ii; = ~ pg 


V 


dxo 


= 0 


La troisieme equation est triviale (0 = 0), on peut done tirer de ce systeme deux 
relations pertinentes : 


f a n "f ^12 — 0 
\a 12 + b 22 = -pg 

La deuxieme equation d'equilibre est celle qui concerne la condition limite sur la 
surface immergee que l'on a ecrite : 

= p e gx 2 ■ T x sur S t0A 

Le premier terme de cette equation vaut: 


(Til 0"i 2 (T 13 


-1 


-<Tn" 

<t 12 °22 °23 


0 

= _ 

<t 12 

. a 13 °23 ^33- 


. 0 . 


. <7 13. 


On peut done ecrire : 

r— (Tn . Pe9 x 2 

] ~ a i 2 - 0 sur S tA 

{ -(Tig = 0 


La troisieme equation est triviale (0 = 0). On remplace les termes de <r par leurs 
expressions dans les deux premieres equations : 


f-UnXi - b ±1 x 2 = p e gx 2 

{ -ai 2 Xi - ^12*2 = 0 


sur S t0A 


Or sur la surface S t0A on peut ecrire x t = 0, et il vient done : 

(~b 1± x 2 = p e gx 2 
l ~b 12 x 2 — 0 

Et finalement, on a les deux relations suivantes : 


[Kx = -Pe9 

l ^12 = 0 
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La derniere equation d'equilibre dont on dispose est celle qui concerne la surface a fair 
libre notee S tB , que l'on a ecrite : 


O 




= 0 sur S t0B 


Le premier terme de cette equation s'ecrit: 


O 




On o 12 

<*13 


' 1/V2 ' 

1 

"On — <*12" 

<*12 <*22 
.<*13 <*23 

<*23 

<*33. 


-1/V2 
. 0 

= vf 

<*12 — <*22 
.<*13 — <*23. 


Par consequent, l'equation d'equilibre se traduit par 

f<*n ~~ <*12 = 0 
Jcri 2 -a 22 = 0 surS t0B 

(<*13 — 0 2 3 = 0 

La troisieme equation de ce systeme est triviale une fois de plus, et les deux autres 
nous permettent d'ecrire : 

fan*! + b ± ix 2 - <*12*1 - ^12*2 = 0 

1<*12^<1 "b b 12 X 2 ~~ a 22 X l ~ b 2 2 X 2 = 0 t0B 

Or sur la surface S t0B on peut ecrire x x — x 2 , et ces termes peuvent done se simplifier 
dans le systeme precedent. 11 vient done : 

( a n ^11 — a 12 ~ b X2 — 0 
(.<*12 "b b x2 — cl 22 — b 22 — 0 

En recapitulant les trois equations d'equilibre, on obtient un systeme avec six equations 
et six inconnues : 


<*n + b 12 — 0 
<*12 + b 22 = -pg 
< b ±1 = —p e g 

b 12 — 0 

<*11 "b ^11 — <*12 — b X2 — 0 
^<*12 "b b 12 — a 22 — b 22 — 0 

En travaillant dans ce systeme pour calculer toutes les variables, et en tenant compte 
du fait que l'on a demontre plus haut que o 33 = v(a 1± + o 22 ), on peut definir 
entierement notre champ de contraintes en precisant les huit coefficients inconnus : 
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( a ±1 = 0 et b tl = -p e g 

I 0 x 2 = “PeP et *>12 = 0 
j a 22 = (p - 2p e )p et b 22 = (p e - p)g 
Va 33 = v(p - 2 p e )p et b 33 = -vpg 

et done : 

Oil = 011*1 + *>11*2 = ~P e P*2 
cji2 = a 12 Xx + b 12 x 2 = -p e gx t 
0X3 = 0 

(722 = a 22 *l + b 22 x 2 = (p - 2 p e )gx ± + (p e - p)px 2 

°23 — 0 

a- 33 = a 33 x x + b 33 x 2 = v(p - 2 p e )gxx - vpgx 2 

Le tenseur de contrainte est done entierement connu. Les trois valeurs non-nulles du 
champ de deformation se deduisent ensuite de la loi de Hooke : 

'\_ 'Y % 

£n — (on ~ V(J 22 ~ ’t' cr 33 ) = — (ctxx — va 22 — va 33 ) + — (b 1:L — vb 22 — vb 33 ) 

1 + v 1 + v 1 + v 

£ 12 = ^ cr 12 = ^ a 12*l ”1 £ bi 2 X 2 

1 

£ 22 = — ( _ 'V<Tii + CT22 — V<T 33 ) 

= y (-VCLXX + a 22 - va 33 ) + -j (-vh xl + b 22 - vb 33 ) 

Tous calculs faits, on obtient: 

—vg g 

£ll = *i ~y~ (1 + V) (p - 2 p e ) + X 2 - (1 + V)(vp - p e ) 

1 + v 

£ 12 “ ~ x lPed £ 

£22 = *11 (1 - v 2 )(p - 2p e ) + x 2 1 [p(v 2 - 1) + p e (l + v)] 

Le champ de deformation est done lui aussi entierement defini. 
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6. La seule contrainte tangentielle non nulle est la contrainte a 12 , et cette contrainte 
augmente (en valeur absolue) proportionnellement avec la coordonnee x 1 . Elle est 
done maximale (en valeur absolue) au point B, ou elle vaut: 

0i2 = -p e gh 
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